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La presente invention concerne un precede et un dispositif de 
simulation interactive du contact entre au moins un premier objet 
5 deformable et au moins un second objet avec un pas predetermine de 
temps d'echantillonnage d'un modele simule. 

On a deja propose d'effectuer une simulation de mesures 
d'interpenetration entre un objet rigide et un objet deformable a partir 
d'estimations volumiques ou de distance, notamment pour des 
10 applications de chirurgie virtuelle ou un outil virtuel rigide de chirurgie 
coopere avec un organe virtuel deformable du corps humain. 

Toutefois, selon ces methodes la relation entre la mesure 
d'interpenetration et les forces de contact de reaction n'ont pas de base 
physique et des forces artificielles peuvent s'appliquer a des nceuds des 
15 maillages des objets qui ne sont pas en contact, ce qui nuit a la fiabilite, 
ou les forces de contact ne repondent pas aux conditions du probleme de 
Signorini. 

La presente invention vise a remedier aux inconvenients precites et 
a permettre de realiser une simulation interactive en temps reel du contact 

20 entre des objets, dont au moins certains sont deformables, de fagon 
simplifiee et economique tout en respectant les contraintes des lois 
physiques qui regissent les contacts, de telle maniere que les contacts 
simules entre des objets soient fiables et qu'ainsi la stabilite de la 
simulation soit garantie. 

25 Ces buts sont atteints grace a un procede de simulation interactive 

du contact entre au moins un premier objet deformable et.au moins un 
second objet avec un pas predetermine de temps d'echantillonnage d r un 
modele simule, 

caracterise en ce que : 

30 (a) on calcule au prealable les parametres decrivant les caracteristiques 
physiques de chacun des objets, teis que la geometrie et la 
mecanique des materiaux de chacun des objets, et on stocke ces 
parametres dans une memoire, 

(b) au debut de chaque pas de temps d'echantillonnage du modele 
35 simule, on procede au niveau de chaque objet a une analyse en 
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temps reel du comportement propre de I'objet pour predire les 
positions, vitesses et accelerations de cet objet selon un mouvement 
libre qui ne tient pas compte d'eventuels contacts ulterieurs, 

(c) a chaque pas de temps d'echantillonnage du modele simule, on 
analyse en temps reel, au niveau d'une scene globale comportant les 
objets susceptibles de venir en contact, des paires d'objets qui sont 
detectes en intersection, et on etablit une liste de groupes de 
collisions qui contient une chaine d'objets en collision et une 
description des collisions, 

(d) a chaque pas de temps d'echantillonnage du modele simule, on 
rapatrie en temps reel, pour chaque groupe de collisions, des 
parametres representant les caracteristiques physiques des objets et 
la description des collisions, de maniere a determiner, pour chaque 
cas, la solution au probleme de Signorini qui gouverne le contact 
entre deux objets dans le cas d'un glissement relatif pur, 

(e) a la fin de chaque pas de temps d'echantillonnage du modele simule, 
on procede au niveau de chaque objet a une visualisation en temps 
reel du comportement propre de I'objet suite a la collision, et 

(f) I'ensemble des traitements en temps reel s'effectue avec un pas de 
temps de calcul plus court que le pas de temps d'echantillonnage du 

modfle-shmiter-deH^ oCl 



I'utilisateur peut intervenir directement en cours de simulation. 
Lors de I'etape a) de calcul prealable des parametres decrivant les 
caracteristiques physiques de chacun des objets, on utilise pour les 
parametres decrivant la mecanique des materiaux une description des 
deformations de type elements finis, avec le remplissage et l'inversion de 
matrices, la resolution de systemes d'equation et le stockage des donnees 
en memoire. 

Selon un mode particulier de realisation, chaque objet est decnt 
dans une configuration au repos comme un ensemble de triangles 
reproduisant sa surface et un ensemble de tetraedres decrivant I'interieur 
de I'objet. 
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Avantageusement, chaque triangle est decrit par trois points, places 
dans un ordre qui permet de calculer des normales qui sont 
invanablement dirigees vers I'exterieur de I'objet. 

De preference, les deformations des objets sont interpolees par la 
5 methode des elements finis en utilisant un maiiiage tetraedrique lineaire. 

A chaque pas de temps de calcul, on integre lors de I'etape b) au 
niveau d'un objet, les forces explicites appliquees sur I'objet, qui sont deja 
connues au depart du pas de calcul, de maniere a definir le mouvement 
qu'elles creent sur I'objet, tandis que la valeur des forces de contact 
10 implicites, qui dependent elles-memes du mouvement des objets dans le 
pas de temps de calcul, est determinee lors de I'etape d) de recherche au 
niveau d'une scene globale, de la solution au probleme de Signorini. 

Lors de I'etape c) d'analyse au niveau d'une scene globale, on detecte 
geometriquement les intersections existantes entre les objets de la scene 
15 pour extraire des couples d'elements d'objets en intersection, une 
longueur et une direction d'interpenetration entre les deux elements d'un 
couple d'elements d'objets. 

Selon une variante de realisation, lors de I'etape c) d'analyse au 
niveau d'une scene globale, pour extraire des couples d'elements d'objets 
en intersection, une longueur et une direction d'interpenetration entre les 
deux elements d'un couple d'elements d'objet, on prend egalement en 
compte un mouvement intermediate des objets entre le pas de calcul 
precedent et le pas de calcul courant, pour calculer une direction 
d'interference privilegiee entre ces objets. 

Lors de I'etape d) de recherche de la solution au probleme de 
Signorini, on reconstruit les points extremes d'application de la force de 
contact entre deux objets soumis a une collision lorsque ces points 
extremes d'application n'ont pas ete determines a I'etape precedente. 

Selon un mode de realisation particulier, lors de I'etape d) lors de 
I'etape d) de recherche de la solution au probleme de Signorini, dans le 
cas d'une intersection segment-segment de deux objets en triangle, les 
deux points choisis pour constituer les points extremes d'application de la 
force de contact entre les deux objets soumis a une collision sont situes a 
('intersection de chacun des deux segments avec le plan forme par la face 
du triangle en intersection. 
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Selon un autre mode de realisation particulier, lors de I'etape d) de 
recherche de ia solution au probleme de Signorini, dans le cas d'une 
intersection point-face de deux objets en triangle, un premier point choisi 
pour constituer un point extreme duplication de la force de contact entre 
5 les deux objets soumis a une collision est le point de Intersection tandis 
que le second point extreme duplication de la force de contact entre les 
deux objets soumis a une collision est la projection du premier point 
extreme sur la face du triangle en intersection. 

Selon un aspect particulier de la presente invention, on utilise les 
10 coordonnees barycentriques pour repartir les deplacements et les forces 
des points ^application de la force de contact entre les points extremes 
d'application de la force de contact en effectuant une interpolation lineaire 
pour une moderation en elements finis. 

De facon plus particuliere, on peut calculer la distance 5 
15 d'interpenetration entre les deux points extremes d'application de la force 
de contact dans le cas d'un contact segment-segment entre un premier 
segment et un second segment d'un second triangle a partir de I'equation 
suivante : 



5 = [a, b ; c,] [[a 1 - «][^] - * 1 



■p] 



(1) 



ou : a et 1-a sont les coordonnees baryceTitriqaes-sur-ie-premier- 

segment, 

p et 1-p sont les coordonnees barycentriques sur le deuxieme 
segment, 

25 a, bi q sont les coordonnees de la direction nj d'interpenetration, 

Wi et W 2 sont les coordonnees du premier segment, 
Vi et V 2 sont les coordonnees du second segment. 
On peut egalement calculer la distance 8 d'interpenetration entre 
les deux points extremes d'application de la force de contact dans le cas 
30 d'un contact point-plan entre un point d'un second triangle et un plan d'un 
premier triangle a partir de I'equation suivante : 
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ou : a, p et y sont les coordonnees barycentriques sur le premier triangle, 

a\ b\ q sont les coordonnees de la direction n ( d'interpenetration, 

5 Wi, W2, W 3 sont les coordonnees du premier triangle, 

Vi represente les coordonnees du point de contact constitue par un 
sommet du second triangle. 

Apres avoir determine les points d f application des forces de contact 
entre deux objets en collision, lors de Petape d) on procede au transfert 
10 des caracteristiques mecaniques des objets dans Pespace defini des 
contacts dans lequei on traite Pensemble d'un groupe de rn contacts avec 
n objets oumetn sont des entiers. 

De facjon plus particuliere, lors de Petape d) on considere la masse 
et Pinertie d'un objet de fagon globale, en son centre de gravite et on 
15 etablit une relation instantanee entre les forces de contact f c dans la 
direction du contact, les accelerations 8^ dues aux contraintes dans la 
meme direction et les accelerations libres 5" libre dans la meme direction 
connues lors de Petape c) au niveau d'une scene globale, selon ('equation 
suivante : 

20 

j c M- i j c T f c +<5'; bro (3) 

ou : J c est une matrice jacobienne m*6n qui transfere le mouvement 
instantane lineaire et angulaire dans Pespace des contacts, 

25 3 T c est la matrice transposee de 3 C , 

M est une matrice bloc diagonale correspondant a la masse et 
Pinertie de n objets du groupe de contacts. 

Lors de Petape d) pour le transport des caracteristiques mecaniques 
locales, on etablit une relation entre : 
30 • la difference de depiacement (Ur) des points du maillage 

deformable representant Pobjet i a Pinstant k, entre la 
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deformation libre (U'idibre) et la deformation contrainte (U^c) soit 

U k = U'|c,c - U'uibre 

. Les positions relatives des objets, libres et contraintes, dans 
l'espace des contacts : Snbre et 8 C . 

6 = £l=i n N c ! U k ! + 8| ib re (4) 

ou N c ' est une matrice de passage de l'espace des deplacements du 
maillage vers l'espace des deplacements aux contacts. 

De meme, on etablit une relation entre les forces dans l'espace des 
contacts f c et les forces dans l'espace des forces de deformation F k . 

F k = (N c ') T / c (5) 

De facon plus particuliere, lors de I'etape d) on etablit une relation 
lineaire instantanee de caracterisation des deformations ou deplacement 
de contact 5 C a partir des forces de contact f c et des deplacements libres 
6 libre dus aux mouvements libres integrant uniquement les forces connues 
de facon explicite en debut de pas de temps de calcul, selon I'equation 
suivante : 

20 S 0 = n; a (u k ,) (N c y ] / c + \ bre (6) 

otn— H 1 — BtnoneTTratrioe-nde-^^ 

maillage vers l'espace des deplacements de contacts, 

(N^.) T est la matrice transposee de N^, 

A est une matrice permettant de definir la deformation de I'objet au 
niveau local, de telle sorte que si U k represente le vecteur du 
deplacement dans le repere local de I'objet a I'instant courant et U k -i 
represente le vecteur de deplacement dans le repere local de I'objet 
au pas de calcul precedent dont les valeurs instantanees sont 
connues au debut du pas de calcul courant, on a : 

U K = A (U k -i) F k + b (U k -i) (7) 

ou F k est un vecteur representant les forces exterieures appliquees a 
I'objet exprimees dans le repere local, et b est un vecteur qui a une 
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valeur dans Pespace des deplacements et qui depend du modele de 
deformation de Pobjet 

Dans un cas plus general, lors de Petape d) on etablit une relation 
instantanee de caracterisation des deformations ou deplacements de 
5 contact 5 C a partir des forces de contact f c et des deplacements libres Sabre 
dus aux mouvements libres integrant uniquement les forces connues de 
fagon explicite en debut de pas de temps de calcul, selon Pequation 
suivante : 

10 S c = [e * 2 j c m- 1 j c T +i: i= ; n; a (u w ) (n;> t ] / 0 + 8 ubre (8) 

ou : J c est une matrice jacobienne m*6n qui transfere le mouvement 
instantane lineaire et angulaire dans Tespace des contacts, 

J I est la matrice transposee de J c f 

M est une matrice bloc diagonale correspondant a la masse et 
Tinertie des n objets du groupe de contacts, 

6 est une constante dependant de la methode d'integration en 
temps, 

N' c est une matrice de passage de Pespace des deplacements du 
maillage vers Pespace des deplacements de contacts, 

(W c ) y est la matrice transposee de N' c , 

A est une matrice permettant de definir la deformation de Pobjet au 
niveau local, de telle sorte que si represente le vecteur du 
deplacement dans le repere local de Pobjet a Pinstant courant et Ur-i 
represente le vecteur de deplacement dans le repere local de Pobjet 
au pas de calcul precedent dont les valeurs instantanees sont 
connues au debut du pas de calcul courant, on a : 

U K - A (U k -i) F k + b (U k -i) (7) 

F k est un vecteur representant les forces exterieures appliquees a 
Pobjet exprimees dans le repere local, et 
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b est un vecteur qui a une valeur dans I'espace des (Replacements et 
qui depend du modele de deformation de I'objet. 
On peut ainsi remarquer que dans ie cas d'un objet rigide 
indeformable, on a U k = U k -i qui traduit I'absence de modification en 
fonction du temps du vecteur U k . 

Avantageusement, le procede selon ('invention comprend en outre 
une etape de couplage avec un module d'interface haptique pour produire 
un retour de sensation haptique sur un dispositif mecanique par iequel un 
operateur manipule les objets dans une scene virtuelle. 

L'invention concerne egalement un dispositif de simulation 
interactive du contact entre au moins un premier objet deformable et au 
moins un second objet avec un pas predetermine de temps 
d'echantilionnage d'un modele simule, caracterise en ce qu'il comprend : 

(a) un module de calcul prealable des parametres decrivant les 
caracteristiques physiques de chacun des objets, tels que la 
geometrie et la mecanique des materiaux de chacun des objets, 

(b) une memoire de stockage des parametres prealablement calcules 
dans le module de calcul, 

(c) un module de couplage avec une interface utilisateur comprenant un 
rjigpogitif mecanique tenu par un utilisateur lui permettant d'exercer 

virtuellement des forces sur lesdits objets dans une scene du modele 
simule, 

(d) un ecran de visualisation pour visualiser lesdits objets represents 
sous forme de maillages, 

(e) une unite centrale de traitement associee a des moyens d'entree, 
comprenant au moins 

el) un module d'analyse d'objet pour analyser en temps reel au 
niveau de chaque objet le comportement propre de I'objet pour 
predire les positions, vitesses et accelerations de cet objet selon 
un mouvement fibre qui ne tient pas compte d'eventuels 
contacts ulterieurs, 

e2) un modele d'analyse d'une scene globale comportant les objets 
susceptibles de venir en contact, pour analyser en temps reel 
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des paires d'objets qui sont detectes en interaction et etablir 
une iiste de groupes de collisions qui contient une chaine 
d'objets en collision et une description des collisions, 

e3) un module de rapatriement en temps reel, pour chaque groupe 
5 de collisions, des parametres representant les caracteristiques 

physiques des objets et la description des collisions pour 
determiner, pour chaque cas, la solution au probleme de 
Signorini qui gouverne le contact entre deux objets dans le cas 
d'un glissement relatif pur, 

10 e4) un module de traitement de chaque objet pour permettre en 

temps reel au niveau de chaque objet une visualisation en 
temps reel du comportement propre de Tobjet suite a une 
collision, et 

e5) des moyens de determination d T un pas de calcul plus court que 
15 le pas de temps d'echantiilonnage du modele simule de maniere 

a definir une simulation interactive. 

Avantageusement, le dispositif comprend des moyens de production 
d r un retour de sensation haptique sur Interface utilisateur. 

Selon une caracteristtque avantageuse, le pas de calcul correspond 
20 a une frequence egale ou superieure a environ 500 Hz. 

D r autres caracteristiques et avantages ressortiront de la description 
suivante de modes particuiiers de ['invention, donnes a titre d'exemples, 
en reference aux dessins annexes, sur lesquels : 

- la Figure 1 est un schema montrant les differentes etapes d'un 
25 procede de simulation interactive du contact entre objets, selon 

['invention, 

- la Figure 2 est un schema montrant differents niveaux de 
traitement de Pinteraction entre objets au cours de differentes etapes du 
procede de simulation de la Figure 1, 

30 - fes Figures 3A a 3C represented trois exemples d'interaction entre 

deux objets virtuels representees par des triangles, 

- la Figure 4 represente une interaction entre deux objets en 
triangle dans le cas d'une intersection segment-segment, 

- la Figure 5 represente une interaction entre deux objets en 
35 triangle dans le cas d'une intersection point-face, 
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- la Figure 6 illustre schematiquement le cas d'une collision entre 
deux objets pour laquelle la configuration peut etre definie a partir de 
seuls criteres geometriques, 




- la Figure 7 illustre schematiquement le cas d'une collision entre 
deux objets pour laquelle la configuration est definie en tenant compte 
d'un mouvement intermediaire, 

- la Figure 8 est un schema-bloc montrant les constituants de base 
d'un dispositif de simulation interactive du contact entre objets, selon 
I'invention, 

- la Figure 9 montre un exemple de contact entre un objet virtuel 
deformable et un autre objet virtuel, et 

- les Figures 10A a 10C montrent trois positions relatives differentes 
entre un objet virtuel deformable, en forme de pince, et un objet virtuel 
rigide lors d'un processus de mise en place de I'objet virtuel deformable en 
forme de pince sur I'objet virtuel rigide. 

La Figure 8 illustre de facon schematique un exemple de dispositif 
permettant de mettre en reuvre I'invention et de realiser la simulation 
interactive en temps reel du contact entre objets tout en permettant 
notamment d'avoir un retour de sensation haptique. 

Une unite centrale de traitement 100, qui peut etre constitute a 
partir d'un ordinateur de type classique, permet d'effectuer les differents 
calculs n ecessaires pour realiser une simulation. 

Un ecran de visualisation 107 relie a I'ordinateur 100 par une 
interface graphique permet 1'affichage d'objets represents sous la forme 
de maillage comprenant des nceuds ou sommets reliant des segments ou 
aretes. 

Des informations peuvent etre fournies a I'ordinateur 100 a partir 
d'une interface utilisateur classique 103 pouvant comprendre un clavier et 
par exemple une souris et constituant des moyens d'entree. 

Un dispositif mecanique specifique 104 tenu par un utilisateur relie 
par un module de couplage 101 a I'ordinateur 100 peut en outre etre 
prevu pour permettre a I'utilisateur d'exercer virtuellement des forces sur 
les objets dans une scene d'un modele simule. Un tel dispositif mecanique 
104 et le module de couplage 101 constituent une interface haptique qui 
permet a I'utilisateur d'exercer une sollicitation sur les objets virtuels de la 
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scene et de recevoir en retour une simulation haptique qui est une 
reponse fournie par la simulation du contact entre objets. 

Uordinateur 100 comprend de fagon ciassiquement au moins un 
processeur, une memoire permanente de stockage de programmes et de 
5 donnees et une memoire de travail cooperant avec le processeur. Des 
supports externes de memoire (disquettes, CD-ROM,...) ou un modem de 
liaison a un reseau peuvent naturellement etre utilises pour charger dans 
Pordinateur des programmes ou donnees permettant d'effectuer tout ou 
partie du traitement de simulation. Sur la Figure 8, on a implement 

10 represents a titre symbolique un exempie de memoire de stockage 102 
cooperant avec le module 100 et pouvant etre constitute par Tun ou 
Pautre des types de memoires sus-mentionnes. 

D T une maniere generate, au debut d'un processus de simulation les 
parametres decrivant la geometrie et la mecanique des materiaux des 

15 objets a simuler sont calcules dans Punite centrale 100 et stockes dans 
une zone memoire de la memoire 102. 

Pour caracteriser les deformations mecaniques des objets, il est 
utilise lors du traitement par Punite centrale 100 une description des 
deformations de type elements finis. Ceia se traduit par le remplissage et 

20 Pinversion de matrices, la resolution de systemes d'equations et ie 
stockage des donnees dans la memoire 102 associee a Punite centrale 
100. 

Les positions et formes courantes des objets sont evaluees en 
fonction des sollicitations exercees et des lois mecaniques qui regissent les 

25 objets de la scene. 

Selon ('invention, pour garantir une simulation stable, il est pris en 
compte dans le calcul des contacts simules entre objets les lois physiques 
qui regissent le contact. Pour permettre une simulation en temps reel, 
c'est-a-dire avec un delai tres court et borne separant la sollicitation 

30 exercee par Putilisateur via Pinterface haptique 104, 101 et la reponse 
fournie par Punite centrale de traitement 100 a cette interface haptique 
104, 101, le dispositif de simulation met en oeuvre trois modules 
principaux qui sont sollicites iterativement a chaque pas de temps 
d'echantillonnage du modele simule. Par ailleurs, Pensemble des 

35 traitements en temps reel s f effectue avec un pas de temps de calcul plus 
court que le pas de temps d'echantillonnage du modele simule. 
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Les trois modules principaux du dispositif de simulation mettant en 
osuvre les differentes etapes du procede de simulation se presentent 
essentiellement de la facon suivante : 

un premier module dit "mecanique" se situant au niveau de chaque 
5 objet et decrivant son comportement propre, permet de faire evoluer 

la position et la forme de I'objet suivant les forces et les iieux des 
forces exerces. Ce module est appele en debut de pas de calcul pour 
predire les positions, vitesses et accelerations des objets sans tenir 
compte du contact puis sera de nouveau mobilise pour tenir compte 
10 des forces calculees dans un troisieme module dit "de traitement du 

contact". 

Un second module dit de "detection de collision", se situant au 
niveau de la scene globale etablit des paires d'objets qui sont 
detectes en intersection. Ce module, de facon optionnelle, peut creer 

15 des mouvements intermediaires entre les pas de la simulation pour 

savoir quand et comment les objets sont entres en intersection. Ce 
module est regi avant tout par des lois geometriques optimisees qui 
permettent d'accelerer le calcul afin d'obtenir une chaine d'objets en 
collision et une description des collisions. Un groupe de collision est 

20 ainsi un ensemble d'objets relies entre eux par au moins une 

collision. Un objet entre dans un groupe s'il est en collision avec au 

moins un des objets du groupe. Une collision est decrite 

obligatoirement par la paire d'objets en collision et par le lieu de la 
collision en utilisant soit les elements geometriques de base (par 

25 exemple deux triangles ou deux surfaces) en intersection, soit par un 

segment reliant les deux points qui sont localement les plus 
interpenetres. 

Un troisieme module dit "traitement du contact" est appele par le 
module de "detection de collision" et fait appel en retour au module 

30 "mecanique". Pour chaque groupe de collisions, le module de 

traitement du contact rapatrie les caracteristiques physiques des 
objets et la description des collisions. Le module est apte a 
determiner, pour chaque cas, la solution au probleme de Signorini 
qui gouverne le contact entre deux objets dans le cas du glissement 

35 pur. 
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^invention permet d'effectuer une simulation interactive. Une 
simulation est definie par ie pas de temps d'echantillonnage du modele 
simule et par le pas de temps de calcul de ce modele. Le precede selon 
Tinvention met en osuvre un pas de temps de calcul qui est toujours 
5 inferieur au pas de temps choisi pour Pechantillonnage, ce qui permet 
d'avoir une simulation interactive ou I'utilisateur va pouvoir intervenir 
directement en cours de simulation. 

La Figure 1 resume les principals etapes du procede selon 
('invention qui met en oeuvre une boucle de simulation utilisant les trois 
10 modules principaux precites instailes dans Pordinateur 100 de la Figure 8. 
La Figure 2 illustre les differents niveaux de traitement entre objets au 
cours des differentes etapes du procede de simulation. 

Une premiere etape de traitement 130 utilise le module dit 
"mecanique" et se situe au niveau de chaque objet (niveau objet 3). Les 
15 informations sont fournies a travers un module de couplage 120 a partir 
de interface utilisateur 110 ou interface haptique qui determine la 
position et la forme de chaque objet (informations 135 elaborees a Petape 
130). 

Dans cette premiere etape 130, une moderation 13 prend en 

20 compte chaque objet ou outil 201, 202, 203 de fagon individuelle sans 
tenir compte d'interactions ulterieures eventuelles et permet de faire 
evoluer la position et la forme de Pobjet suivant les forces et les lieux des 
forces exerces a partir de Pinterface utilisateur 110. 

Une deuxieme etape du traitement 140 utilise le module dit 

25 "detection de collision" et se situe au niveau d r une scene globale (niveau 
scene 4). Les informations 135 elaborees a Petape 130 sont utilisees a 
Petape 140 pour etabiir des paires d T objets qui sont detectes en 
intersection. Au cours de cette etape 130, il est elabore une liste de 
groupes de collisions (informations 145) contenant une chafne d f objets en 

30 collision et une description des collisions. Dans cette etape 140, une 
modeiisation 14 prend ainsi en compte une paire d'objets en intersection 
tels que les objets 201, 202 au niveau d r une scene globale, 

Une troisieme etape de traitement 150 utilise le module dit 
"traitement du contact" et se situe au niveau d ! une scene globale (niveau 

35 scene 5). Les informations 145 elaborees a Tetape 140 ainsi que les 
informations 135 elaborees a Tetape 130 sont utilisees pour determiner 
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pour chaque cas la solution au probleme de Signorini qui gouveme le 
contact entre deux objets dans le cas de glissement pur (information 155). 
Dans cette etape 150, une modelisation 15 prend ainsi en compte 
I'interaction entre deux objets tels que les objets 201, 202 au niveau d'une 
5 scene globale, les caracteristiques physiques des objets et la description 
des collisions etant rapatriees pour chaque groupe de collisions par le 
module "traitement de contact". 

La troisieme etape de traitement 150 fournit des informations 155 
concernant des forces et des lieux qui sont transmises au premier module 
10 "mecanique" au cours d'une quatrieme etape de traitement 160 qui se 
situe a nouveau au niveau des objets (niveau objet 6). Dans cette etape 
160, le resultat du traitement de simulation en temps reel peut etre 
simplement normalise dans une etape 170 de visualisation ou peut etre 
transmis en retour a travers le couplage 120 vers ['interface utilisateur 110 
15 pour conferer a 1'utilisateur un retour de sensation haptique. Dans cette 
etape finale, une modelisation 16 prend ainsi a nouveau en compte 
chaque objet ou outil 201, 202, 203 individuellement tout en ayant tenu 
compte des contacts precedemment simules. 

Au niveau de I'objet, dans une forme de realisation preferee, on 
20 decrit I'objet comme un ensemble de triangles reproduisant sa surface et 
un ensemble de tetraedres pour decrire son interieur, le tout dans une 
c onfig uration de repos. Cette configuration correspond a la forme de 
I'objet lorsqu'aucune force ne lui est appliquee. 

Avantageusement, les triangles sont decrits par trois points, places 
dans un ordre qui permet au calcul des normales d'etre invariablement 
dirigees vers I'exterieur de I'objet. Les surfaces des objets sont fermees de 
fagon a pouvoir distinguer un exterieur d'un interieur. 

Les deformations des objets sont interpolees par la methode des 
elements finis en utilisant un maillage tetraedrique lineaire. Le dispositif 
permet de simuler differentes lois de comportement sous reserve que I'on 
puisse en extraire localement, approximativement et pour un pas de 
calcul, une relation lineaire entre les forces exercees et les deplacements 
autour d'une configuration locale. 

Si U k represente le vecteur de deplacement dans le repere local 
d'un objet a I'instant courant t et si U k - X represente le vecteur de 
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deplacement dans le repere local de I'objet au pas de calcul precedent M, 
on a la relation suivante : 

U k = A(U W ) F k + b(U M ) (7) 

5 

ou : A est une matrice permettant de definir la deformation de I'objet au 
niveau local, 

Ft< est un vecteur representant les forces exterieures appliquees a 
I'objet exprimees dans le repere local, 
10 b est un vecteur qui a une valeur dans I'espace des deplacernents et 

qui depend du modele de deformation de I'objet, et 
Ur-i est un vecteur dont les valeurs instantanees sont connues au 
debut du pas de calcul, de ('instant courant t. 

15 Avantageusement, une distinction est faite entre le mouvement 

global de I'objet decrit par la relation fondamentale de la dynamique et la 
deformation de I'objet autour d'une configuration courante decrite par la 
loi de deformation, 

Les forces appliquees sur I'objet se distinguent par leur caractere 

20 explicite ou implicite. Une force explicite est d'ores et deja connue au 
depart du pas de calcul, et il faut I'integrer pour connaitre le mouvement 
qu'elle cree sur I'objet Les forces de contact, au contraire, sont irnplicites 
au sens ou elles dependent elles-memes du mouvement des objets dans 
le pas de temps. Pour un pas de temps, on integre done les forces 

25 explicites et on passe alors au niveau de la scene globale pour pouvoir 
trouver la valeur des forces irnplicites. 

Apres ('integration des forces explicites au niveau de chaque objet, 
on appelle position et forme "libra" la configuration des objets dans la 
scene. Cette configuration est obtenue sans ('intervention des forces de 

30 contact. Ainsi, un "mouvement libre" est un mouvement integrant 
uniquement les forces connues de fagon explicite en debut de pas de 
temps de calcul. Par suite on considere ici que le mouvement libre est le 
mouvement obtenu sur un pas de temps quand les forces de contact ne 
sont pas integrees. 

35 Le system e propose met alors en ceuvre un processus de detection 

de collision qui permet de tester geometriquement les intersections 
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existantes entre les objets 223, 230 de la scene et les directions 
privilegiees pour sortir les objets de cette collision (Figures 6 et 7). 

Si I'on considere un objet 221 qui, dans une position 223, vient en 
interaction avec un autre objet 230, la direction privilegiee peut etre 
calculee uniquement sur des criteres geometriques (cas de la Figure 6) ou 
peut dependre de la configuration par laquelle les objets 223, 230 sont 
entres en collision en tenant compte d'un mouvement intermediate 222 
d'au moins Tun des objets entre le pas de calcul precedent et le pas de 
calcul courant (cas de la Figure 7). 

Dans tous les cas le processus de detection des collisions permet 
d'extraire des couples d'elements d'objets en intersection, une longueur et 
une direction d'interference entre ces deux elements. 

Dans le cas privilegie d'une description de la surface des objets par 
des triangles, un element est soit un point, soit un segment, soit la face 
d'un triangle. La detection de collision peut prendre en compte trois cas 
canoniques d'intersection entre deux objets : intersection point/triangle, 
intersection segment/segment, intersection triangle/point. 

Le procede peut egalement fournir un ensemble d'elements d'objet 
en proximite qui pourrait potentiellement entrer en collision suite a 
Integration des forces de contact. Une distance et une direction separant 
ces elements sont alors catculees. 

GrScp. a la descri ption de toutes les interferences et proximites 

entre les objets, I'ensemble des groupes de collisions de la scene peut etre 
construit. Chaque groupe va alors passer dans le module de contact. 

Apres la detection d'une collision, un traitement dans le module de 
contact permet de determiner la configuration du premier contact entre 
deux pas de temps de calcul. 

Si Ton prend en compte une description de la surface des objets par 
des triangles lorsqu'une zone d'interpenetration entre deux objets a ete 
definie a un instant T du pas d'echantillonnage de la simulation, il est 
extrait une liste de triangles composant la paire d'objets. Si Ton a un objet 
rigide et un objet deformable, les coordonnees du triangle representant 
I'objet deformable sont traduites dans le cadre de reference de I'objet 
rigide a des instants distincts d'echantillonnage de simulation T etT-1. 

Pour toute paire possible triangle/triangle il est effectue une 
interpolation lineaire du deplacement de trois points Di, D 2 , D 3 du triangle 
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deformable entre ies positions initiate et finale aux instants discrets 
d'echantillonnage T et T-l. On peut ensuite operer trois types de tests 
differents representee sur Ies Figures 3A a 3C sur la paire comprenant un 
triangle deformable 20 defini par des sommets Di, D 2 , D 3 et presentant 
5 diverses positions successives 21, 22, 23 entre Ies instants T-l et T et un 
triangle rigide 30 defini par des sommets Ri, R 2 , R3. 

Le test 1 (Figure 3A) correspond au cas ou un plan de collision est 
forme par le triangle rigide et va etablir une contrainte sur le point 
concerne du triangle deformable* 
10 Le test 2 (Figure 3B) correspond au cas ou un plan de collision est 

forme par un segment rigide et un segment deformable a Tinstant de 
collision t, et va etablir une contrainte sur deux points de I'objet 
deformable. 

Le test 3 (Figure 3C) correspond au cas ou un plan de collision est 
15 forme par I'objet deformable a I'instant de collision et va etablir une 
contrainte sur trois points du triangle deformable. 

^invention peut s'appliquer aussi bien au cas de collisions entre un 
objet rigide et un objet deformable qu'au cas de collisions entre deux 
objets deformables. 

20 Le module de contact va etre decrit ci~dessous de fagon plus 

detaillee en reference au cas prefere d'une description des objets en 
triangle. 

Le module de contact est appele autant de fois qu'il y a de groupes 
de collisions dans la scene a I'instant du calcul. Le module de detection de 
25 collisions a stocke dans un espace memoire pour chaque collision : 

- la normale, 

- la paire d r objets et Ies elements concernes par la collision, 

- (eventuellement) les points duplication de la force de contact 

Si Talgorithme de detection de collision ne donne pas Ies points 
30 duplication de la force de contact (cas illustre sur la Figure 6 d'une 
detection sans mouvement intermediate), il faut Ies reconstruire et dans 
tous les cas, il faut interpoler ces points duplications par rapport au 
modele de deformation choish 

Dans le cas prefere d'une description des objets en triangle, on 
35 separe le probleme en deux cas : soit on a deux elements segments, soit 
on a un element nceud et un element face. 
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10 



15 



Dans le cas d'une intersection segment/segment entre un triangle 
40 defini par des sommets Pi, P 2 , Ps et un triangle 50 defini par des 
sommets Q t , Q 2 , Qs (Figure 4), les deux points choisis 41, 51 sont situes a 
I'intersection de chacun des deux segments Qi, Q 2 , resp. Pi, P 2 avec le 
plan forme par la face de I'autre triangle en intersection. 

Le vecteur reliant les deux points trouves 41, 51 est appele 5. 
Dans le cas d'une intersection Point/Face entre un triangle 60 defini 
par des sommets Pi, P 2 , P 3 et un triangle 70 defini par des sommets Qi, 
Q 2 , Qs (Figure 5), le premier point choisi est le point de I'intersection et le 
second est la projection de ce point sur la face de I'intersection. Le vecteur 
reliant les deux points trouves est appele 5. 

On peut utiliser I'algorithme d'intersection avec mouvement 
intermediate propose par Xavier Provot (Collision and Self-collision 
handling in cloth model dedicated to design garment. Graphics Interface 
1997, 177-189) qui permet d'obtenir une configuration approchee entre 
les deux triangles au moment de la collision. 

Pour repartir les deplacements et les forces de ces points, dans le 
cas privilegie d'une interpolation lineaire pour la modelisation en elements 
finis (triangles - tetraedres), on utilise les coordonnees barycentriques. 

* Pour pouvoir calculer correctement les deformations sur les objets, 
il faut garantir la non-interpenetration. L'interpolation des elements est 
done utilisee de facon a n'avoir qu'une inconnue de force et de distance 
-p^r-coTTtectr-Pour-si^ 
interpolation lineaire pour les elements finis. 

On a ainsi dans le cas d'un contact segment/segment (Figure 4) la 

relation suivante : 



u alpha et 1-alpha sont les coordonnees barycentriques sur le premier 
segment Qi, Q 2 et Beta, 1-Beta sur le deuxieme segment Pi, P 2 , 
as, bi, q sont les coordonnees de la normale n s du triangle 40, 
Wi, W 2 sont les coordonnees du premier segment Qi, Q 2 , 
Vi, V 2 sont les coordonnees du deuxieme segment Pi, P 2 . 
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Dans le cas cTun contact point/plan (Figure 5), a {/ b\ , q sont les 
coordonnees de la normale n,- du triangle 60, la distance d'interpenetration 
8 entre les triangles 60 et 70 s'ecrit : 
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ou : alpha, beta, gamma sont les coordonnees barycentriques sur ie 
triangle 60 (addition = 1). 
10 Wi, W 2 , W 3 sont les coordonnees du premier triangle 60, et Vi sont 

les coordonnees du point de contact 71 constitue par un sommet Qi 
du second triangle 70, 

On utilise une interpolation identique pour la force de contact. 
Une fois trouve le point duplication des forces de contact, les 
15 caracteristiques mecaniques des objets sont transferees dans Tespace 
defini des contacts. Pour la suite, on suppose que Ton traite un groupe de 
m contacts avec n objets. 

Pour le transport des caracteristiques mecaniques globales dans le 
cas ou Ton considere la masse et Tinertie de I'objet de fagon globale, en 
20 son centre de gravite, on peut utiliser une matrice jacobienne classique, 
definie dans les travaux de Ruspini (Diego Ruspini & Oussama "A 
Framework for Multi-Contact Multi-Body Dynamic Simulation and Haptic 
Display", Proceedings of the 2000 IEEE/RS3 International Conference on 
Intelligent Robots and Systems"). Grace a cette matrice, on a une relation 
25 instantanee entre les forces de contact f c dans la direction du contact, les 
accelerations 8" c (contraintes) et 8'Wes dans la meme direction : 

<5" c = j c M- l j c T f c+ 5 r ; bre (3) 

30 ou J c est une matrice jacobienne m*6n qui transfere le mouvement 
instantane iineaire et angulaire dans Tespace des contacts, M est une 
matrice bloc diagonale correspondant a la masse et Pinertie des n objets 
du groupe de contacts et J C T est la matrice transposee de 3 C . 
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Pour le transport des caracteristiques locales, on peut utiliser 
Interpolation definie pour les triangles. On a ainsi une relation entre les 
deplacements des points du maillage deformable et les deplacements dans 
I'espace des contacts et la meme interpolation peut etre utilisee entre les 
5 forces de contact et les forces sur les noeuds. 

En reprenant la relation lineaire instantanee de caracterisation des 
deformations, on obtient : 

S c = [S,; Nj A (U k ,) (N C ') T ] f c + % bre (6) 
ou Nc 1 represente la matrice de passage de I'espace des deplacements du 
10 maillage vers I'espace des deplacements aux contacts et A est une matrice 
com me defini plus haut. 

La moderation du contact est choisie de facon a respecter au mieux 
les lois physiques. Par souci de simplification on peut neanmoins de 
preference considerer que les contacts ne changeront pas de direction 
15 durant la resolution du calcul meme si en pratique ce n'est pas strictement 
le cas. 

Le premier postulat du probleme de Signorini est qu'il n'y a pas 
d'interference entre les objets s'ils sont solides (les matieres ne se 
melangent pas). Ainsi, on souhaite qu'apres la resolution du probleme le 
20 deplacement au contact soit positif ou nul : 

^>-0— (-9-) 

Le second postulat est que Ton est dans le cas d'un contact sans 
25 frottement, done la force de contact est dirigee selon la normale : 

/>0 (10) 

Le troisieme postulat est que la force de contact est non nulle {f c * 0) si et 
seulement si il y a reellement un contact (<5 C =0). Cela cree une relation 
30 complementaire entre les deux vecteurs : 

8 C ±/ C (ID 

Pour pouvoir resoudre le probleme de Signorini, on utilise le 
35 transport des caracteristiques mecaniques. Les effets des caracteristiques 
locales et globales sont additionnes en integrant les accelerations pour se 
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ramener uniquement a une relation entre emplacement et force dans 
Pespace des contacts. 

Avantageusement, on utilise une methode numerique qui a 
tendance a conserver Penergie comme la methode des trapezes (aussi 
5 nommee Tustin), pour integrer ies accelerations et qui peut s ! exprimer de 
la fagon suivante si Ton considere une grandeur X aux instants t (Xt) et 
t+1 (X t+ i) : 

x t+I =x t +?/2dt (x; + x t+ ;) 
10 x tH ^x; + V2dt(x;'+xj>) (12) 

En utiiisant cette methode numerique et Ies equations (4) et (5), on 
obtient la relation suivante : 

15 S c = [>/4dt 2 j c ]vrV^ (13) 

Le coefficient 1/4 peut etre different si Ton utilise une autre 
methode d'integration des accelerations. 

Si le modele mecanique choisi ne comporte pas de caracteristiques 
20 globales et si la masse et rinertie sont integrees au niveau local dans le 
modele de deformation, on utilise uniquement Pequation (5) qui comporte 
deja implicitement une integration numerique en temps. 

5 C - IX; Nj A (U k ,) (N;) T ] f^ ubre (6) 

25 

Les postulats definis dans le probleme de Signorini et Pequation 
lineaire instantanee creant une relation lineaire entre forces de contact et 
deplacements dans Pespace des contacts permettent de formuier le 
contact sous forme de probleme de complementarite lineaire (LCP). Pour 

30 ce type de problemes, il existe de nombreux algorithmes de resolution 
(voir par exemple Murty, K.G., Linear Complementarity, Linear and 
Nonlinear Programming. Internet Edition 1997) qui sont capables de 
resoudre le probleme en un temps compatible avec Ies performances 
demandees par Phaptique. A titre d'exemple, on peut effectuer un calcul a 

35 une frequence de Pordre de 500 Hz a 1000 Hz pour un nombre 
raisonnable de contacts (par exemple 30 a 40 contacts) en utiiisant 
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I'algorithme du Pivot principal sur un ordinateur de type PC Pentium IV 
2GHz. 

La Figure 3 illustre ("interaction entre un objet deformabie 80, tel 

qu'une pince avec un autre objet 90. Dans ce cas, I'objet deformabie 80 
5 est virtuellement attache a Interface haptique (des lors qu'il est tenu par 

I'utilisateur) dans une zone d'un noeud O definissant un repere Ox 0 yoZo. 

Pour de tels points on peut appliquer les conditions de Dirichlet 

A chaque pas de simulation, la configuration de mouvement libre 

donne un espace de contact, Une resolution LCP donne les forces de 
10 contact f { non nulies et un mouvement contraint en est deduit. Ces forces 

(illustrees par la normale U, de coordonnees ai h q sur la Figure 9) sont 

transportees au point O pour creer la force et le couple sur I'interface 

haptique. 

Les Figures 10A a 10C montrent I'exemple d'un objet deformabie 80 
15 constitue par un clip mis en place sur un tube 91. On voit differentes 
deformations du clip 80 dans des positions 81, 82, 83 differentes par 
rapport au tube 91. 
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REVENDICATIONS 

1. Precede de simulation interactive du contact entre au moins un 
5 premier objet deformable (201 ; 20 ; 40 ; 60 ; 80) et au moins un second 
objet (202 ; 30 ; 50 ; 70 ; 90 ; 91) avec un pas predetermine de temps 
d'echantillonnage d'un modele simule, 

caracterise en ce que : 

(a) on calcule au prealable les parametres decrivant ies caracteristiques 
10 physiques de chacun des objets / tels que la geometrie et la 

mecanique des materiaux de chacun des objets, et on stocke ces 
parametres dans une memoire, 

(b) au debut de chaque pas de temps d'echantillonnage du modele 
simule, on procede au niveau de chaque objet a une analyse en 

15 temps reel du comportement propre de l'objet pour predire ies 

positions, vitesses et accelerations de cet objet seion un mouvement 
libra qui ne tient pas compte d'eventuels contacts ulterieurs, 

(c) a chaque pas de temps d'echantillonnage du modele simule, on 
analyse en temps reel, au niveau d'une scene globale comportant les 

20 objets susceptibies de venir en contact, des paires d'objets qui sent 

detectes en intersection, et on etablit une liste de groupes de 
collisions qui contient une chaine d'objets en collision et une 
description des collisions, 

(d) a chaque pas de temps d'echantillonnage du modele simule, on 
25 rapatrie en temps reel, pour chaque groupe de collisions, des 

parametres representant les caracteristiques physiques des objets et 
la description des collisions, de maniere a determiner, pour chaque 
cas, la solution au probleme de Signorini qui gouverne le contact 
entre deux objets dans le cas d r un glissement relatif pur, 

30 (e) a la fin de chaque pas de temps d'echantillonnage du modele simule, 
on procede au niveau de chaque objet a une visualisation en temps 
reel du comportement propre de Tobjet suite a la collision, et 
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(f) Tensembie des traitements en temps reel s'effectue avec un pas de 
temps de calcul plus court que le pas de temps d'echantillonnage du 
modele simule, de maniere a definir une simulation interactive ou 
I'utilisateur peut intervenir directement en cours de simulation. 

5 2. Procede selon la revendication 1, caracterise en ce que lors de 

I'etape a) de calcul prealable des parametres decrivant les caracteristiques 
physiques de chacun des objets, on utilise pour les parametres decrivant 
la mecanique des materiaux une description des deformations de type 
elements finis, avec le remplissage et I'inversion de matrices, la resolution 

10 de systemes d'equation et le stockage des donnees en memoire. 

3. Procede selon la revendication 1 ou 2, caracterise en ce que 
chaque objet est decrit dans une configuration au repos comme un 
ensemble de triangles reproduisant sa surface et un ensemble de 
tetraedres decrivant I'interieur de I'objet. 

15 4. Procede selon la revendication 3, caracterise en ce que chaque 

triangle est decrit par trois points, places dans un ordre qui permet de 
calculer des normaies qui sont invariablement dirigees vers I'exterieur de 
I'objet 

5. Procede selon la revendication 3 ou 4, caracterise en ce que les 
20 deformations des objets sont interpolees par la methode des elements 
" finis en utilisant un maiitalpTeWeTfr^ ; 

6. Procede selon I'une quelconque des revendications 1 a 5, 
caracterise en ce que a chaque pas de temps de calcul, on integre lors de 
I'etape b) au niveau d'un objet, les forces explicites appliquees sur I'objet, 

25 qui sont deja connues au depart du pas de calcul, de maniere a definir le 
mouvement qu'elles creent sur I'objet, tandis que la valeur des forces de 
contact implicites, qui dependent elles-memes du mouvement des objets 
dans le pas de temps de calcul, est determinee lors de I'etape d) de 
recherche au niveau d'une scene globale, de la solution au probleme de 

30 Signorini. 

7. Procede selon I'une quelconque des revendications 1 a 6, 
caracterise en ce que lors de I'etape c) d'analyse au niveau d'une scene 
globale, on detecte geometriquement les intersections existantes entre les 
objets de la scene pour extraire des couples d'elements d'objets en 
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intersection, une longueur et une direction d'interpenetration entre les 
deux elements d'un couple d'elements d'objets. 

8. Procede selon la revendication 7, caracterise en ce que lors de 
Petape c) d'analyse au niveau d'une scene globale, pour extraire des 
5 couples d'elements d'objets en intersection, une longueur et une direction 
d'interpenetration entre les deux elements d'un couple d'elements d'objet, 
on prend egalement en compte un mouvement intermediate des objets 
entre le pas de calcul precedent et le pas de calcul courant, pour calculer 
une direction d'interference privilegiee entre ces objets. 

10 9. Procede selon Tune quelconque des revendications 3 a 5 et selon 

la revendication 7, caracterise en ce que lors de Petape d) de recherche de 
la solution au problerne de Signorini, on reconstruct les points extremes 
duplication de la force de contact entre deux objets soumis a une 
collision lorsque ces points extremes duplication n'ont pas ete 

15 determines a Petape precedente. 

10. Procede selon la revendication 9, caracterise en ce que lors de 
Petape d) de recherche de la solution au problerne de Signorini, dans le 
cas d f une intersection segment-segment de deux objets en triangle (40, 
50), les deux points (41, 51) choisis pour constituer les points extremes 

20 duplication de la force de contact entre les deux objets (40, 50) soumis a 
une collision sont situes a ('intersection de chacun des deux segments 
(P1P2, Q1Q2) avec le plan forme par la face du triangle en intersection. 

11. Procede selon la revendication 9, caracterise en ce que lors de 
Petape d) de recherche de la solution au problerne de Signorini, dans le 

25 cas d f une intersection point-face de deux objets en triangle (60, 70) , un 
premier point (71) choisi pour constituer un point extreme duplication de 
la force de contact entre les deux objets (60, 70) soumis a une collision 
est le point de ['intersection tandis que le second point extreme 
duplication de la force de contact entre les deux objets soumis a une 

30 collision est la projection (61) du premier point extreme (71) sur la face du 
triangle en intersection (60). 

12. Procede selon Pune quelconque des revendications 9 a 11, 
caracterise en ce qu'on utilise les coordonnees barycentriques pour 
repartir les displacements et les forces des points duplication de la force 
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de contact entre !es points extremes d'application de la force de contact 
en effectuant une interpolation lineaire pour une moderation en elements 
finis. 

13. Procede selon les revendications 10 et 12, caracterise en ce que 
Ton calcule la distance 5 d'interpenetration entre les deux points extremes 
(41, 51) d'application de la force de contact dans le cas d'un contact 
segment-segment entre un premier segment (Qi Q 2 ) et un second 
segment (Pi P 2 ) d'un second triangle a partir de I'equation suivante : 



ou : a et 1-a sont les coordonnees barycentriques sur le premier 
segment (Qi Q 2 ), 

p et 1-p sont les coordonnees barycentriques sur le deuxieme 
segment (Pi P 2 ), 

aj bj q sont les coordonnees de la direction n, d'interpenetration, 
Wi et W 2 sont les coordonnees du premier segment Qi Q 2 , 

V! et V 2 sont les coordonnees du second segment Pi P 2 . 

14. Procede selon les revendications 11 et 12, caracterise en ce que 
Ton calcule la distance 5 d'interpenetration entre les deux points extremes 
(61, 71) d'application de la force de contact dans le cas d'un contact point- 
plan entre un point (71) d'un second triangle et un plan (Pi P 2 P 3 ) d'un 
premier triangle a partir de I'equation suivante : 



6 = [a s bj cj [a I - a] 




- CP 



Ell 



d =[a t b, cj [a jS y] W 2 - V, (2) 

KJ 



a, p et y sont les coordonnees barycentriques sur le premier triangle, 
aj b, q sont les coordonnees de la direction n { d'interpenetration, 
Wi, W 2 , W 3 sont les coordonnees du premier triangle (Pi P 2 P 3 ), 
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Vi represente les coordonnees du point de contact constitue par un 
sommet (Qi) du second triangle (Qi Q 2 Qs). 

15. Precede seion Pune quelconque des revendications 1 a 14, 
caracterise en ce qu'apres avoir determine les points duplication des 

5 forces de contact entre deux objets en collision, lors de Petape d) on 
procede au transfert des caracteristiques mecaniques des objets dans 
Pespace defmi des contacts dans lequel on traite Pensemble d'un groupe 
de m contacts avec n objets oumetn sont des entiers. 

16. Procede selon la revendication 15, caracterise en ce que lors de 
10 Petape d) on considere la masse et Pinertie d T un objet de fagon globale, en 

son centre de gravite et on etablit une relation instantanee entre les forces 
de contact f c dans la direction du contact, les accelerations 8^ dues aux 
contraintes dans la meme direction et les accelerations libres 8), te dans la 
meme direction connues lors de Petape c) au niveau d'une scene gfobale, 
15 selon ('equation suivante ; 

J C M-Vf c+ S'; bra ( 3) 

ou : J c est une matrice jacobienne m*6n qui transfere le mouvement 
20 instantane lineaire et angulaire dans Pespace des contacts, 

J I est la matrice transposee de J c , 

M est une matrice bloc diagonale correspondant a la masse et 
Pinertie de n objets du groupe de contacts. 

17. Procede selon Pune quelconque des revendications 9 a 14 et 
25 Pune quelconque des revendications 15 et 16, caracterise en ce que lors 

de Petape d) pour le transport des caracteristiques mecaniques locales, on 
etablit une relation entre : 

® la difference de deplacement (Uk) des points du maillage 
deformable representant Pobjet i a Pinstant k, entre la 
30 deformation libre (UVjibre) et la deformation contrainte (lA,c) soit 
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• les positions relatives des objets, libres et contraintes, dans 
I'espace des contacts : 8| ib re et 8 C , avec 

5 = Xl=l" Nc' U k ' + Siibre (4) 

ou N c j est une matrice de passage de I'espace des deplacements du 
maillage vers I'espace des deplacements aux contacts, et 

on etablit une relation entre les forces dans I'espace des contacts f c 
et les forces dans I'espace des forces de deformation F k :: 

Ftc = (Nc*) t /c (5) 

18. Procede selon la revendication 1, caracterise en ce que lors de 
I'etape d) on etablit une relation lineaire instantanee de caracterisation des 
deformations ou deplacement de contact 8 C a partir des forces de contact 
f c et des deplacements libres 8 iib re dus aux mouvements libres integrant 
uniquement les forces connues de facon explicite en debut de pas de 
temps de calcul, selon I'equation suivante : 

<5 C = [s w n n; a (uj (n;) t ] f c + $ libre (6) 

ou : N' e est une matrice de passage de I'espace des deplacements du 
maillage vers I'espace des deplacements de contacts, 

(.^^)Lest-la-matrice-transposee_de_W4, 

A est une matrice permettant de definir la deformation de I'objet au 
niveau local, de telle sorte que si U k represente le vecteur du 
deplacement dans le repere local de I'objet a I'instant courant et U k -i 
represente le vecteur de deplacement dans le repere local de I'objet 
au pas de calcul precedent dont les valeurs instantanees sont 
connues au debut du pas de calcul courant, on a : 

U K = A (Uk-i) F k + b (U k -i) (7) 

ou F k est un vecteur representant les forces exterieures appliquees a 
I'objet exprimees dans le repere local, et 

b est un vecteur qui a une valeur dans I'espace des deplacements 
et qui depend du modele de deformation de I'objet. 
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19. Precede selon la revendication 17, caracterise en ce que lors de 
Petape d) on etablit une relation instantanee de caracterisation des 
deformations ou deplacements de contact 8 C a partir des forces de contact 
f c et des deplacements libres Snbre dus aux mouvements libres integrant 
5 uniquement les forces connues de fagon explicite en debut de pas de 
temps de calcul, selon ('equation suivante : 

<5 C - [e dr j c M" 1 jj+^j n; a ou (N L y ] / e + 8 Hbre m 

10 ou : J c est une matrice jacobienne rn*6n qui transfere le mouvement 
instantane lineaire et angulaire dans Pespace des contacts, 

est la matrice transposee de 3 C ' 

M est une matrice bloc diagonale correspondant a la masse et 
Pinertie des n objets du groupe de contacts, 

15 0 est une constante dependant de la methode d'integration en 

temps, 

N[, est une matrice de passage de Pespace des deplacements du 
maillage vers Pespace des deplacements de contacts, 

(N<) T est la matrice transposee de IM^, 

20 A est une matrice permettant de definir la deformation de Pobjet au 

niveau local, de telle sorte que si U k represente le vecteur du 
deplacement dans le repere local de Pobjet a Pinstant courant et Uk-i 
represente le vecteur de deplacement dans le repere local de Pobjet 
au pas de calcul precedent dont les valeurs instantanees sont 

25 connues au debut du pas de calcul courant, on a : 

U K - A (U k -i) F k + b (Uk-i) (7) 

ou F k est un vecteur representant les forces exterieures appliquees a 
30 Pobjet exprimees dans le repere local, et 

b est un vecteur qui a une valeur dans Pespace des deplacements 
et qui depend du modele de deformation de Pobjet 
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20. Procede selon I'une quelconque des revendications 1 a 19, 
caracterise en ce qu'il comprend en outre une etape de eouplage avec un 
module d'interface haptique pour produire un retour de sensation haptique 
sur un dispositif mecanique par lequel un operateur manipule les objets 

5 dans une scene virtuelle. 

21. Dispositif de simulation interactive du contact entre au moins 
un premier objet deformable (201 ; 20 ; 40 ; 60 ; 80) et au moins un 
second objet (202 : 30 ; 50 ; 70 ; 90 ; 91) avec un pas predetermine de 
temps d'echantillonnage d'un modele simule, caracterise en ce qu'il 

10 comprend : 

(a) un module (100) de calcul prealable des parametres decrivant les 
caracteristiques physiques de chacun des objets, tels que la 
geometrie et la mecanique des materiaux de chacun des objets, 

(b) une memoire (102) de stockage des parametres prealablement 
15 calcules dans le module (100) de calcul, 

(c) un module- (101) de eouplage avec une interface utilisateur (104) 
comprenant un dispositif mecanique tenu par un utilisateur lui 
permettant d'exercer virtuellement des forces sur lesdits objets dans 
une scene du modele simule, 

20 (d) un ecran de visualisation (107) pour visualiser lesdits objets 
representes sous forme de maillages, 
(e) une unite centrale de traitement (100) associee a des moyens 
d'entree (103), comprenant au moins 

el) un module d'analyse d'objet pour analyser en temps reel au 
25 niveau de chaque objet le comportement propre de I'objet pour 

predire les positions, vitesses et accelerations de cet objet selon 
un mouvement libre qui ne tient pas compte d'eventuels 
contacts ulterieurs, 

e2) un modele d'analyse d'une scene globale comportant les objets 
30 susceptibles de venir en contact, pour analyser en temps reel 

des paires d'objets qui sont detectes en interaction et etablir 
une liste de groupes de collisions qui contient une chaine 
d'objets en collision et une description des collisions, 
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e3) un module de rapatriement en temps reel, pour chaque groupe 
de collisions, des parametres representant les caracteristiques 
physiques des objets et la description des collisions pour 
determiner, pour chaque cas, la solution au probleme de 
5 Signorini qui gouverne le contact entre deux objets dans le cas 

d'un glissement relatif pur, 

e4) un module de traitement de chaque objet pour permettre en 
temps reel au niveau de chaque objet une visualisation en 
temps reel du comportement propre de Pobjet suite a une 
10 collision, et 

e5) des moyens de determination d f un pas de calcul plus court que 
le pas de temps d'echantillonnage du modele simule de maniere 
a definir une simulation interactive, 

22, Dispositif selon la revendication 21, caracterise en ce qu'il 
15 comprend des moyens de production d'un retour de sensation haptique 

sur ('interface utilisateur (104). 

23. Dispositif selon la revendication 21 ou 22, caracterise en ce 
que le pas de calcul correspond a une frequence egale ou superieure a 
environ 500 Hz. 
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